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I. INTRODUCTION 

Symmetry breaking is an ubiquitous feature of the 
low temperature behavior in condensed matter physics. 
Solids or Neel antiferromagnets are phases that break 
some essential symmetries of the physical laws: trans- 
lational symmetry or rotational spin symmetry. Under- 
standing the nature of the broken symmetries, discrete 
or continuous, allows to understand the nature of the 
elementary excitations and to predict the low-energy be- 
havior of the materials (Goldstone modes, Mermin Wag- 
ner theorem, topological defects, ...). In some phases, at 
first glance, the symmetry content may be hidden: as for 
example in Helium liquids. The first obvious character 
is the absence of translation symmetry breaking and ab- 
sence of a solid phase at zero temperature. It was very 
early understood (F. London) that this absence of so- 
lidification is due to the many-body quantum dynamics 
and the Helium phases have been named quantum liq- 
uids to be contrasted to the more "classical liquids". It 
was only decades after discovery of the ''He superfluidity 
that the nature of the order parameter was unveiled. The 
understanding of the '^He order parameter has also been 
heavily dependent on group symmetry considerations. 

A parallel can be developed between this distinction of 
quantum liquids versus classical solids and that of spin 
liquids (SL's) versus Neel ordered phases. Neel ordered 
phases at least break translational symmetry of the lat- 
tice and rotational symmetry of the spins. They can 
be described by a local order parameter and a Landau 
theory, whereas SL do not break any lattice symmetries 
nor spin rotation symmetry and cannot be described by 
a local order parameter. Similarly to ^He, SL's can be 
characterized by an internal hidden, more or less complex 
order. 

In this paper we are mainly concerned with topolog- 
ical SL. These SL's are characterized at T = by ex- 



ponentially decaying correlations for all local observables 
(spins, dimers or spin nematic operators) and a spin gap 
to bulk excitations. They contrast to critical SL which 
have algebraic correlations and gapless excitations. It 
has been understood very early ' that the elementary 
excitations of these resonating valence bond (RVB) SL 
carry a spin-^ contrarily to the spin-1 magnons of the 
Neel antiferromagnets. These emergent excitations are 
called spinous. A natural framework to describe the SL 
physics is the use of effective theories with the fractional 
particles as elementary building blocks (parton construc- 
tion). Going from the original spins to these fractional- 
ized spinous implies the introduction of gauge fields in 
which the spinous are deconfined (SL) or glued (Neel or- 
der). At first glance these approaches introduce via the 
gauge fields a considerable (infinite) amount of degrees 
of freedom. In fact the number of possible distinct SL's 
is limited by the requirement that their physical observ- 
ables do not break any lattice or spin symmetry and the 
enumeration of the different classes of distinct SL's can 
be done through group theory analysis. 

This was understood ten years ago by X.-G. Wen who 
developed a classification of symmetric SL using Projec- 
tive Symmetry Group technique (PSG). ' The analysis of 
Wen for fermionic spinous on the square lattice was ex- 
tended by Wang and Vishwanath to bosonic spinous. In 
these works, the definition of a SL is limited to spin sys- 
tems that do not break any symmetry, neither SU (2) spin 
symmetry nor lattice symmetries nor time reversal sym- 
metry. These SL's have been dubbed by Wen symmetric 
SL. This definition excludes chiral SL which break time- 
reversal symmetry (and some minimal amount of lattice 
symmetry) but which do not break SU(2) and do not 
have long range order in spin-spin or dimer-dimer corre- 
lations. 

In the wake of Laughlin theory of FQHE, chiral 
SL's have been very popular at the end of the eighties 
(Kalmeyer and Laughlin Wen et al. ' , Yang et al. ), 
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but, in the absence of indisputable candidates, this op- 
tion has nearly disappeared from many discussions in the 
last decade. 

Non-planar structures are quite ubiquitous in classi- 
cal frustrated magnetism,'' and are associated to scalar 
chirality: Si ■ {S2 x S3) ^ 0. In some cases where the 
ground state is non-planar this chirality can persists at 
finite temperature ' although the magnetic order itself 
is absent for T > (Mermin- Wagner) . A similar phe- 
nomenon may take place in quantum systems at T = 0. 
There, the usual scenario is that of a gradual reduction 
of the Neel order parameter when the strength of the 
quantum fluctuations is increased. At some point the 
sublattice magnetization vanishes and the SU (2) sym- 
metry is restored (leading to a SL). Now, if the ordered 
magnetic structure is chiral, the time-reversal symmetry 
T may still be broken at the point where the magnetic 
order disappears, hence leading to a time-reversal sym- 
metry breaking (TRSB) SL. Some TRSB SL have in- 
deed been recently proposed on the kagome lattice^^'^'^ 
and there are probably other examples.'^''"' 

The goal of this paper is to revisit the PSG analysis by 
relaxing the time-reversal symmetry constraint in order 
to include chiral SL's. The framework used here is the 
Schwinger-boson mean-field theory (SBMFT). But, as 
for the symmetric PSG, the symmetry considerations we 
use here should also be valid to classify SL in presence of 
moderate fluctuations beyond mean-field. 

The paper is organized as follow. Sections II and III 
are reviews, to keep this article self-contained. Section II 
is a description of SBMFT to fix the notations and pre- 
cise the present understanding of this approach. Sec. Ill 
starts by recalling the gauge invariance of SBMFT and 
then describes how the PSG is used to enforce the SL's 
symmetries on mean-field theories. 

In Sec. IV, the concept of PSG is extended to include 
all chiral SL's. In Sec. V all the chiral and non chiral 
SL theories with explicit nearest neighbor gauge fields 
on the triangular lattice are derived. As an example of 
application we propose a chiral SL as the ground state 
of a ring-exchange model on the triangular lattice. The 
physical meaning of the fiuxes and their expressions in 
terms of spin operators is developed in Sec. VI, as well as 
the question of topological loops on finite size samples. 
Sec. VII is the conclusion. Appendices contain proofs of 
some statements in the main text, technical details and 
further applications to the square and kagome lattices. 
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II. SCHWINGER BOSON MEAN-FIELD 
THEORY (SBMFT) 

We consider a spin Hamiltonian -ffo({Si}i=i...7Vs) on 
a periodic lattice with Ng spins, each of length S. Hq 
can contain Heisenberg interaction or more complicated 
terms such as cyclic exchange, all invariant under global 
spin rotations {SU(2) symmetry) and by time-reversal 
transformation T (i/o({S,}) = i/o({-SJ)). We insist 
on these symmetries since they are the basis of our con- 
struction. 

Finding the ground state (GS) of a quantum spin prob- 
lem is notoriously difhcult problem and the SBMFT pro- 
vides an approximate way to treat the problem. This ap- 
proach can be summarized by the following steps: i) The 
spin operators (hence the Hamiltonian) are expressed us- 
ing Schinwer bosons, ii) A suitable rotationally-invariant 
mean-field decoupling leads to a quadratic Hamiltonian 
Hmf- iii) ^^MF is diagonalized using a Bogoliubov trans- 
formation and solved self-consistently. 

A. Bosonic operators and bond operators 

Let m to be the number of sites per unit cell in the 
lattice, and Nm number of unit-cells, so that Ng — 
is the total number of sites. We define the two bosonic 
operators that create a spin a = ±1/2 (or cr =t or |) 
on site i. The spin operators read: 

St=%h, (lb) 

s- = bi;b,^. (ic) 

The Hamiltonian is thus a polynomial of bosonic opera- 
tors with only even degree terms. These relations imply 
that the commutation relations [5f,5f] — ie^^^Sf are 

verified. As for the total spin, it reads Si — ^ -1- l), 

where rii = bLbt^ + bl^hi is the total number of bosons 
at site i. To fix the "length" of the spins, the following 
constraint must therefore be imposed on physical states: 

= = (2) 

(7 

In traditional MF theories, the MF parameter is the 
order parameter (as for example the magnetization (8,^)) 
and the MF Hamiltonian consequently breaks the initial 
Hamiltonian symmetries, except in the high temperature 
phase where the MF parameter is zero. Here, we would 
like to describe SL's that do not break any symmetry. 
Thus we are going to express Hq using quadratic bosonic 
operators, requiring their invariance by global spin rota- 
tions. 

The expectation value of these operators will then be 
used as mean-field parameters, insuring that the MF 



Hamiltonian respects the rotational invariance. Only lin- 
ear combinations of the two following operators and of 
their hermitian conjugates obey this property: 

^ij = ^(Mii-Mjt)' (3a) 

B^, = libltbn + blM)- (3b) 

i and j are lattice sites and these operators are thus bond 
operators. They are linked by the relation 

: BjjB,, : +AI^A,, = ^n,{n, - S,,) (4) 

where : . : means normal ordering. 

Any Hamiltonian invariant by global spin rotation can 
be expressed in terms of these operators only. For ex- 
ample, an Heisenberg term Si ■ Sj where i ^ j can be 
decoupled as 

Si ■ Sj — : BjjBij : —Al^Aij, (5a) 
= 2 : %B,, : (5b) 
= S^~2Al^A,. (5c) 

where the first line is true whatever the boson number, 
but the last two lines use Eq. 4 and suppose that the 
constraint of Eq. 2 is strictly respected. 

To make clear the physical significance of these two 
bond operators in the case S = 5, we write them in 
terms of projection operators Pg on the singlet state and 
Pt on the triplet states: 

44 = (6a) 
■■Bj^B,, : = liPt^Ps). (6b) 

We see in Eq. 6, that : BjjBij : represents a ferromagnetic 

contribution to Eq. 5a, whereas A\jAij gives the singlet 
contribution. 

B. The mean- field approximation 

We now need two successive approximations to obtain 
a quadratic and solvable Hamiltonian. We first relax the 
constraint on the boson number by imposing it only on 
average: 

{n.) = n. (7) 

where n does not need to be a integer. To implement this 
constraint, a Lagrange multiplier (or chemical potential) 
\i is introduced at each site i and the term \i{n — 
rii) is added to the Hamiltonian. k can be continuously 
varied to interpolate between the classical limit (k = 00) 
and the extreme quantum limit (k — > 0). 
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It should be reminded in general that fixing k = 25* to 
study a spin-5 model is not necessarily the best choice as 
in the SBMFT (Sf) = |k(k + 2).1'' An alternative choice 
could be to fix k, in such a way that the spin fluctuations 
and not the spin length have the correct value. 

In a second step, bond operators fluctuations are ne- 
glected and a MF Hamiltonian -//mf that is linear in bond 
operators is obtained. For instance: 



{Bl^)B,,+Bl^{B,,)~m,)\^. (8) 



We replace (Bij) and {Aij) by complex bond parameters 
Aij and Bij. This MF approximation can be seen as the 
first term of a large N expansion of a Sp{N) theory. ' ' 
The steps are explained in details in Ref. .im in the very 
similar case of an SU{N) theory. This zero'th order 
expansion can be pursued to the first order. The MF 
Hamiltonian is now a quadratic bosonic operator. It can 
be written in terms of a 2Ns x 2Ns complex matrix M 
and of a real number cq depending on the Aij and Bij 
and on the Lagrange multipliers A^: 



where 



tt 



MF 



£0- 



\ If 2t' 



(9) 



The ex- 



pression for M and eo depend on Hq and on the chosen 
decoupling (for example using Eq. 5a, 5b or 5c). 

The set of mean- field parameters {Aij, Bij} appearing 
in Hmf is called an Ansatz. Up to an equivalence rela- 
tion that will be described in the next section, an Ansatz 
defines a specific phase (ground state and excitations). 
Depending on the value of k, this state can either have 
Neel long range order, or the bosons are gapped (several 
types of SL are then possible). 

In the following we will explain and exploit the relation 
which exists between regular classical magnetic orders ' 
and SL's. 

To enforce self-consistency, the following conditions 
should be obeyed: 



Aij {-^ij 

which are equivalent to 



and Bij = 



B, 



dFMF „ J QFmf 
and 



^A^ 



dBi, 



0, 



(10) 



(11) 



where Fmf is the MF free energy, together with the con- 
straint 



(12) 



The next step is to calculate the mean values of the op- 
erators Aij and Bij either in the GS of HyiF if the tem- 
perature is zero, or in the equilibrium state for non-zero 
temperatures. In both cases one needs to use a Bogoli- 
ubov transformation to diagonalize Hmf- As this trans- 
formation is often explained in the simple case of 2 x 2 
matrices (or for particular sparse matrices), we explain 
the algorithm in a completely general case in App. A. 



C. Choice of bond fields: Aij and Bij or Aij or Bij 
only. 

As in the example of Eq. 5, the relation 4 can be used to 
eliminate Aij or Bij from Hmf ■ If we choose to keep only 
the Bij parameters, M is block diagonal with two blocks 
of size Ns and the vacuum of bosons is a GS. To obey the 
constraint of Eq. 2, we have to adjust the boson densities 
by filling some zero-energy mode(s), therefore breaking 
the SO (3) symmetry. The GS is thus completely clas- 
sical. On the contrary, we can keep the Aij only, but 
then the singlet weight is overestimated, which can be 
introduce some bias on frustrated lattices where short- 
distance correlations are not coUinear. Keeping Aij only 
is a widespread practice in the litterature, but Trumper 
et al.' ' have explicitly shown that the bandwidth of the 
spectrum of excitations of the Heisenberg model on the 
triangular lattice is twice too large when using Aij fields 
only. On the other hand, the use of both Aij and Bij 
restores the correct bandwidth and a improves quanti- 
tatively the excitation spectrum. Note that even on the 
square lattice the simultaneous use of both bond opera- 
tors improves the ground state energy.'" 

From a different point of view Flint and Coleman" 
advise the use of both fields in order to have a large- 
N limit where spin generators are odd under the time- 
reversal symmetry, as it is the case for SU{2). 



III. THE SEARCH OF SL 

Even when considering an Hamiltonian with nearest 
neighbor interactions only, the dimension of the MF pa- 
rameters manifold is exponentially large. Moreover the 
Lagrange multipliers Xi make the search of the station- 
ary points of the MF free energy difficult (constrained 
optimization) as for each considered Ansatz, all Xi must 
be adjusted to calculate the MF free energy. In Ref. 22 
this optimization was carried out (without any simpli- 
fying/symmetry assumption) on square and triangular 
lattices with up to 36 sites. In almost all cases the 
MF ground-state turned out to be highly symmetric, as 
expected, but excited mean-field solutions are however 
highly inhomogeneous (and often not understood yet). 
The problem can be considerably simplified if we restrict 
our search to states respecting some (or all) the symme- 
tries of Hq . Such symmetries are divided into global spin 
rotations, lattice symmetries and time reversal symme- 
try. We have assumed from the beginning that Hq is 
invariant by global spin rotations and chosen the MF ap- 
proximation in such a way that it remains true for Hmf, 
but the choice of a specific Ansatz may or may not break 
other discrete symmetries. The fore-coming section ex- 
plains how to find all Ansatze such as the physical quan- 
tities are invariant by all the lattice symmetries X, either 
strictly (for symmetric SL's) or only up to a time-reversal 
transformation (Chiral SL's). 
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We will now define some groups specific to an Ansatz: 
the invariance gauge group in Sec. Ill A and the projec- 
tive symmetry group in Sec. IIIB. Then, in Sec. IIIC, we 
define the algebraic projective symmetry group, which is 
associated to a lattice symmetry group and not specific 
to a particular Ansatz on this lattice. 

A. Gauge invariance, fluxes and invariance gauge 
goup (IGG) 



Eq. 13 with 9{i) = ir for all lattice sites i. In the par- 
ticular cases where we can divide the lattice in two sub- 
lattices such as Aij — whenever i and j are in the 
same sublattice (bipartite problem), the IGG is enlarged 
to U{1). The later situation corresponds, for instance, 
to an Ansatz on a square lattice with only first-neighbor 
Aij. The transformations of the IGG are then given by 
d{i) — 9 on a. one sublattice and 6{i) — —9 on the other, 
with arbitrary 9 G [0, 27r[. 



Let Q ^ U{1)^'' be the set of gauge transformations. A 
gauge transformation is characterized by an angle 9{i) € 
[0, 2n[ at each site and the operator G which implements 
the associated gauge transformation 

b,,^%,e''"^3'^=G%M (13) 

is given by 

G = cxp ^z^&]J,„(?(j)j . (14) 

A wave function respects a symmetry F if all the 
physical observables measured in the state F\4)) are iden- 
tical to those measured in \4>). It does not mean that 
|(/)) — F\(l))^ but that the two wave functions are equal up 
to a gauge transformation: 3G € Q ,\4)) — GF\4)). 
The action of G on the Ansatz is: 

\S,fc^6,,e'(-'^«+«W), ^^^^ 

such as HyiF remains unaffected by G. We note that 
{Ajk) and {Bjk) are gauge dependent: they are not phys- 
ical quantities as they do not preserve the on-site boson 
number. As any such quantity, their mean values cal- 
culated using H(j is zero when the average is taken on 
all gauge choices. Using Hmf, it can be non-zero as the 
gauge symmetry is explicitly broken by the choice of the 
Ansatz. 

We have seen that changing the gauge modifies the 
Ansatz but not the physical quantities. Conversely, if two 
MF Hamiltonians give rise to the same physical quanti- 
ties, then their Ansatze are linked by a gauge transforma- 
tion. In fact two types of physical quantities are directly 
related to the Ansatz: the MF parameter moduli (related 
to the scalar product of two spins), and the fluxes. The 
fluxes are defined as the arguments of Wilson loop oper- 
ators such as {B.ijBjkBki) or of {A\jAjkA\iAu) . By con- 
struction these quantities are gauge invariant and define 
the Ansatz up to gauge transformations. The physical 
meaning of fluxes will be addressed in Sec. VI. 

The gauge transformations that do not modify a spe- 
cific Ansatz form a subgroup of Q called the invariance 
gauge group (IGG). It always contains the minimal group 
Z2 formed by the identity and by the transformation 



B. The projective symmetry group (PSG) 

Let X be the group of the lattice symmetries of the 
Hamiltonian Hq (translations, rotations, reflections...). 
From now on, for the sake of simplicity, we discard the 
hat on the gauge and symmetry operators. The effect of 
an element X oi X on the bosonic operators is 

X -.b,, ^hx(^)„. (16) 
The effect of X on the Ansatz is: 

{ Ajk Ax(j)X(k), ryj^ 
I ^3k Bx{j)X{k)- 

We know that a gauge transformation does not change 
any physical quantities. What about the lattice sym- 
metries ? We know from Sec. Ill A that if the Ansatze 
before and after the action of X have the same physical 
quantities, they are linked by a gauge transformation: it 
thus exists at least one gauge transformation Gx such 
as GxX leaves the Ansatz unchanged. The set of such 
transformations of Q x X is called the projective symme- 
try group (PSG) of this Ansatz. Note that this group 
only depends on the Ansatz and on A", but not on the 
details of the Hamiltonian. Thus, an Ansatz is said to 
respect a lattice symmetry X if it exists a transformation 
Gx e G such that the Ansatz is invariant by GxX. 

The IGG of an Ansatz is the PSG subgroup formed 
by the set of gauge transformations G/ associated to the 
identity transformation / of A". For each lattice symme- 
try X € X respected by the Ansatz, the set of gauge 
transformations Gx such as GxX is in the PSG is iso- 
morph to the IGG: for any G/ in the IGG, {GiGx)X is 
in the PSG. Thus, the condition for an Ansatz to respect 
all the lattice symmetries is that its PSG is isomorph to 
IGGxX. 



C. The algebraic projective symmetry groups 

An Ansatz is characterised (partially) by its IGG and 
its PSG. In turn, we know from these groups which lat- 
tice symmetries it preserves. Reversely we now want to 
impose lattice symmetries and find all Ansatze that pre- 
serves them. To reach this goal, we proceed in two steps. 
The first one is to find the set of the so-called algebraic 
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PSG's. '' ' They are subgroups oi Q x X verifying alge- 
braic conditions necessarily obeyed by a PSG. Contrary 
to the PSG of an Ansatz, the algebraic PSG's exist inde- 
pendently of any Ansatz and only depend on the lattice 
symmetry group X and on the choice of an IGG (chosen 
as the more general). An algebraic PSG does not depend 
on the details of the lattice such as the positions of the 
sites. However, depending on these details, an algebraic 
PSG may have zero, one, or many compatible Ansatze. 
The second step consists, for a given lattice, in finding 
all the Ansatze compatible with a given algebraic PSG. 

Let us detail the algebraic conditions verified by the 
algebraic PSG's. The group X is characterized by its 
generators xi. . .Xp. A generator Xa has an order Ua G N* 
such as x"" is the identity (if no such integer exists, we 
set Ua = oo). For any transformation X € X, there 
exists a unique ordered product X = with 
< ka < ria ii fia is finite, /ca S ^ if not. The rules used 
to transform an unordered product into an ordered one 
are the algebraic relations of the group. Each of these 
rules implies a constraint on the Gx^ (chosen as one of 
the gauge transformation associated to Xa)- Basically, it 
states that if a lattice symmetry X can be written in sev- 
eral ways using the generators, the gauge transformation 
Gx is independent of the writing (up to an IGG trans- 
formation). The subgroups oi Q x X respecting all these 
constraints are the algebraic PSG's. 

To illustrate the idea, let us consider a basic exam- 
ple where X is generated by two translations xi and 
X2. Both transformations have an infinite order ni = 
77.2 = oo. We have X & X written as product of gen- 
erators X = j;^'ia;™"a;"='a;™''... and we would like o 
write it as X = x\^x^. The need algebraic relation 
is simply the commutation between the two translations 
: 2:1X2 — X2X1. We then have pi = mi -I- 7713 -I- . . . 
and p2 = rn2 + rui -I- . . . . We will now see that this 
implies a constraint on G^^ and G^a • Suppose that we 
have an Ansatz unchanged by Gx^^xi and Gx2X2- Then 
the inverses x^^G~^ or X2^G~^ too are in the PSG. So, 
the product GxiXiGx2X2Xi^G~lx2^G~l G PSG. This 
product has been chosen to make the algebraic relation 
xia;2 = X2X1 (<^ a;iX2a;^^X2^^ = /) appear after the fol- 
lowing manipulations: 

GxiXiGx2X2Xi ^Gx^x2 ^Gx2 € PSG 
^ Gx,{xiGx2Xi^)xiX2X^^X2Hx2G-^^X2^)G-^ e PSG 
^ GxdxiGx2x^'){x2Gx^'x2')G-l e PSG. 

The expressions in parenthesis in the last line are pure 
gauge transformations and the full resulting expression is 
a product of gauge transformations. Thus, we can more 
precisely write: 

GxAxiGx2X^'){x2Gx^X2')G-^ e IGG. (18) 

If the IGG is Z2, this constraint can be written in term 
of the phases 9x(i) of the gauge transformation Gx as: 

dx, (i) + 0x2 {xi^i) - Ox, {x2^i) - 0x2 (0 = P^, (19) 



with p = or 1. This constraint coming from the com- 
mutation relation between xi and X2 must be obeyed by 
all algebraic PSG's. 

It is useless to list all algebraic PSG's for the sim- 
ple reason that some of them are equivalent and give 
Ansatze with the same physical observables. Two (al- 
gebraic or not) PSG's are equivalent if they are related 
by a gauge transformation G: for any gauge transforma- 
tion Gx associated to the lattice symmetry X in the first 
PSG, GGxG~^ belongs to the set of gauge transforma- 
tions associated to X in the second PSG. We are only 
interested in equivalence classes of PSG's. 

Taking algebraic PSG's in different classes does not 
imply that they have no common Ansatze: a trivial ex- 
ample is the Ansatz with only zero parameters, belonging 
to any algebraic PSG's. But each class includes Ansatze 
that are in no other class and have specific physical prop- 
erties. 

Once all the algebraic PSG's classes are determined, it 
remains to find the possible compatible Ansatze for one 
representant of each class. As an example of compati- 
bility condition, let's take the case where X belongs to 
the considered algebraic PSG {i.e. Gx — I). Then an 
Ansatz can be compatible with this algebraic PSG only 
if, for any couple of sites (?,j), Aij — Ax{i)x(j)- If such 
compatible Ansatze exist, they respect the lattice sym- 
metries by construction (in the sense that their physical 
quantities do so). We now want to impose the time rever- 
sal symmetry: among the compatible Ansatze, we only 
keep those that are equivalent to a real Ansatz up to a 
gauge transformation. We call them strictly symmetric 
Ansatze {weakly symmetric ones are defined in the next 
section). 

To completely define an Ansatz, it is sufficient to give 
the algebraic PSG and the values of the MP parameters 
on non symmetry-equivalent bonds. For example, on a 
square (or triangular or kagome) lattice with all usual 
symmetries (see Fig. 2) and only first neighbor interac- 
tions, the Aij and Bij of one bond are enough. 



IV. FROM CHIRAL LONG RANGE ORDERS 
TO CHIRAL SL'S 

We will now show that the zoo of Neel LRO obtained 
from the strictly symmetric Ansatze misses the chiral 
states which are exact ground states of a large num- 
ber of frustrated classical models. This will lead us in 
a straightforward manner to the construction of chiral 
algebraic PSG's in which time reversal and some lattice 
symmetries can be broken (Sec. IV B). This generalised 
framework will then be illustrated on the triangular lat- 
tice in Sec. V and on the square and kagome lattice in 
App. D. 
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A. SU(2) symmetry breaking of symmetric Ansazte 

To simplify, we suppose that all lattice sites are equiv- 
alent by symmetry and only consider Ansatze such as the 
\i are all equal to a single A. Even if an Ansatz is strictly 
symmetric, it does not always represent a SL phase. As 
is well known in SBMFT, a Bose condensation of zero en- 
ergy spinons can occur and leads to Neel order. We will 
discuss how the Ansazte symmetry constraints the mag- 
netic order obtained after condensation, and establish a 
relation with the regular states introduced in Ref. . 

The Bogoliubov bosons creation operators are linear 
combinations of the bi^ and h\^, such as their vacuum 
|0) is a GS of -ffiviF (see App. A). If the GS is unique, it 
must respect all the Hamiltonian symmetries and conse- 
quently, cannot break the global spin rotation invariance. 
But when k increases (we continuously adapt the Ansatz 
to K so that the self-consistency conditions remains ver- 
ified and the PSG remains the same), some eigen en- 
ergie(s) decrease(s) to zero. The GS is then no more 
unique as the zero mode(s) can be more or less popu- 
lated and the phases of each zero mode are free. It is 
then possible to develop a long range spin order. 

This phenomena occurs when no A verifies condition 
12. If A increases the mean number of boson per site 
increases up to a maximal number Kmax- At this point, 
some eigen energies become zero. Increasing A further is 
not possible as the Bogoliubov transformation becomes 
unrealizable (the M matrix of Eq. 9 has non-positive 
eigenvalues). To reach the required number of boson 
per site, we have to fill the zero energy modes b\, b\, 

. . . using coherent states e"i''i+"2''2+- - 16) for example. In 
the thermodynamical limit the fraction of missing bosons 
is macroscopic and a Bose condensation occurs in each 
of the soft modes. The choice of the weight ai of these 
modes fixes the direction of the on-site magnetization. 
Detailed examples of magnetization calculations in a con- 
densate are given by Sachdev.- ' 

In the classical limit (k — >■ oo), all bosons are in the 
condensate and contribute to the on-site magnetization 
mi. The modulus |mi| should be equal to k/2 to satisfy 
Eq. 9. The bi„ operators acquire a non-zero expectation 
value {bi„) and are (up to a gauge transformation) linked 
to mi by : 

\ ) - \ ^«^e'A.g(™r+.™n ) , (20) 

where Arg is the argument of the complex number and 
m^'^'^ are the magnetization components. These values 
are constrained by the Ansatz through: 

- li{blt)ib,t) + {bl){b,l))- (21b) 
The supplementary constraint reads: 

|m,| ~ k/2 (22) 



This extra constraint can make the classical limit prob- 
lem unsolvable: no classical magnetization pattern is 
then compatible with the Ansatz. An example of such 
a situation was studied by Wang and Vishwanath ' (see 
App. C). 

We can take the problem of the classical limit from the 
other side. We begin from a classical state, from which 
we calculate {bia) and the Ansatz (using Eq. 20 and 21). 
What are the conditions on the classical state for the as- 
sociated Ansatz to be strictly symmetric ? As we look 
for an Ansatz respecting all lattice symmetries, the ro- 
tationally invariant quantities (as the spin-spin correla- 
tions) must be invariant by all lattice symmetries, what 
severely limits the classical magnetization pattern. Such 
a state is called a S'0(3)-regular state. Mathematically, 
a state is said to be S'0(3)-regular if for any lattice sym- 
metry X there is a global spin rotation Sx € 50(3) such 
as the state is invariant by SxX. Moreover, the time 
reversal symmetry (i.e. the Ansatz can be chosen to be 
real) imposes the co-planarity of the spins. The set of 
coplanar S'0(3)-regular states can be sent on the set of 
condensed states of strictly symmetric Ansatze. In the 
same way, we define the 0(3)-regular states by including 
global spin flips Si — )> —Si in the group of spin transfor- 
mations. These 0(3)-regular states are listed in Ref. for 
several two-dimensional lattices. The 0(3)-regular states 
are divided in coplanar S'0(3)-regular states and in chiral 
states. In a chiral state, the global inversion — )■ —Si 
cannot be "undone" by a global spin rotation. Equiv- 
alently, there exist three sites i, j, k such as the scalar 
chirality Si ■ (Sj A S^) is non zero: the spins are not copla- 
nar. Then a strictly symmetric Ansatz, upon condensa- 
tion, can only give coplanar S'0(3)-regular states in the 
classical limit, therefore missing all chiral 0(3)-regular 
states. 

This limitation can seem unimportant as most of the 
usual long-range ordered spin models have planar GS's. 
But some new counter examples have recently been 
discovered. The first example is the cyclic exchange 
model on the triangular lattice^" with a four sublattice 
tetrahedral chiral GS (see Fig. 1). More recently, two 
twelve sublattice chiral GS's, with the spins oriented to- 
wards the corners of a cuboctahedron, were discovered 
on the kagome lattice with first and second neighbor 
exchanges^^'" (studied in App. D 2). A systematic study 
of the classical GS's of simple models on different lattices 
has indeed revealed that the GS's are chiral for large 
ranges of interaction values. ' 

The theory of symmetric PSG is unable to encompass 
such chiral states. In the following subsection, we will 
build TRSB SL Ansatze which include, upon condensa- 
tion, all classical regular chiral states. This method was 
already applied to the kagome lattice with up to third 
neighbor interactions, leading to the surprising result of 
a chiral state even in the purely first neighbor model. 
If this state is physically relevant or not is still an open 
question, but independently, it shows that the omission 
of chiral Ansatze has prevented the discovery of more 
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FIG. 2: (Color online) Generators of the lattice symmetries 
X on the triangular and square lattices. Vi is a translation, 
CT is a reflection and TZi is a rotation of order i. 

competitive MF solutions. 



B. The chiral algebraic PSG's: how to include 
weakly symmetric states 

The time-reversal transformation T acts on an Ansatz 
by complex conjugation of the MF parameters. ' If an 
Ansatz respects this symmetry, it is sent to itself by T 
(up to a gauge transformation). So, in an appropriate 
gauge, all parameters can then be chosen real. In most 
previous SBMFT studies, the hypothesis of time reversal 
invariance of the GS was implicit, as only real Ansatze 
were considered. In contrast to SU{2) global spin sym- 
metry that can easily be broken through the Bose conden- 
sation process, no transition is known to produce a chiral 
ordered state out of a T-symmetric Ansatz. Indeed, chi- 
ral Ansatze have loops with complex-valued fluxes which 
evolve continuously with n. We do not expect any sin- 
gular behavior of these (local) fluxes when crossing the 
condensation point, so the generic situation is that a chi- 
ral LRO phase will give rise to a TRSB SL'^-'^ ' when de- 
creasing K. It is of course possible that the lowest-energy 
Ansatz changes with k but such a first-order transition 
has not reason to coincide with the onset of magnetic 
LRO. 

To obtain all chiral SL's we have to explicitly break 
time-reversal symmetry at the MF level, in the Ansatz. 
For 50(3) classical regular states, a lattice transforma- 
tion from X is compensated by a global spin rotation 
(that leaves the Ansatz unchanged). For 0(3) classical 
regular states, a lattice transformation X d X is com- 
pensated by a global spin rotation possibly followed by 
an inversion — >■ — S^. This defines a parity ex to be 
if no spin inversion is needed, and —1 otherwise. In 
a chiral SL, the parity will be deduced from the effect 
of X on the fluxes: ex = 1 if they are unchanged, — 1 
if they are reversed. With this distinction in mind we 
will call weakly symmetric Ansatze (WS) the Ansatze re- 
specting the lattice symmetries up to T, whereas the the 
Ansatze respecting strictly all lattice symmetries and T 
have already been called strictly symmetric (SS) Ansatze 
(all lattice symmetries are even). 

The distinction between even and odd lattice symme- 
tries (as defined by ex) is the basis of the construction 
of all WS Ansatze via the chiral algebraic PSG's. Let us 



consider Xg the subgroup of transformations of X that 
can only be even. Mathematically, X^ is the subgroup 
of X which elements are sent to the identity by all mor- 
phisms from X to Z2 . X^ contains at least all the squares 
of the elements of X as ex^ — e\ = I. But, depending on 
the algebraic relations of A", it may contain more trans- 
formations as we show in the triangular case in Sec. IV C. 
Once Xf. is known, we define the chiral algebraic PSG's of 
X as the algebraic PSG's of X^- The method described 
previously to find all algebraic PSG's applies the same 
way. We define Xo as the set of transformations which 
may be odd {X — Xg^). It contains transformations of 
undetermined parities. 

To filter the weakly symmetric Ansatze from those 
compatible with the chiral algebraic PSG's, we have to 
take care of the transformations of Xq. This gives two 
types of extra constraints. First, same type {A or B) MF 
parameters on bonds linked by such transformation must 
have the same modulus. The second constraint concerns 
their phases, through the fluxes. The phases are gauge 
dependent, but the fluxes are gauge independent. Fluxes 
are sent to their opposite by T and as well as by the 
odd transformations of X . Then are unchanged by even 
transformations. To find all WS Ansatze we then have 
to determine a maximal set of independent elementary 
fluxes and distinguish all possible cases of parities for 
the transformations of Xq {ex = ±1)- 

We can now apply these theoretical considerations to 
find all WS Ansatze on some usual lattices as the tri- 
angular, honeycomb, kagome and square lattice. The 
calculations are detailed for the triangular lattice in the 
following subsections and some results for the kagome 
and square lattice are given in App. D. 



C. Chiral algebraic PSG's of lattices with a 
triangular Bravais lattice 

The first step is to find all chiral algebraic PSG's. As 
already mentioned, they only depend on the symmetries 
of Xe and on the IGG. We choose the most general case 
of IGG~ Z2 and suppose that respects all the lat- 
tice symmetries with the generators described in Fig. 2. 
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These symmetries are those of a triangular lattice, but 
the actual (spin) lattice of Hq can be any lattice with a 
triangular Bravais lattice such as a honeycomb, a kagome 
or more complex lattices. The coordinates (x, y) of a 
point are given in the basis of the translation vectors Vi , 
V2 and the effect of the generators on the coordinates are 



Vi : (x,y) 




f23a) 


V2 : (x,y) 


— !• fx, W + l), 


(23b) 


lZf\ : fx, 


— >■ fx — y, x), 


f23c) 


(7 ' ( X n) 


— >■ f?/ x) 


(23d) 


The algebraic relations in 


A" are: 




V1V2 


= V2V1, 


(24a) 




= / 


(24b) 




- 


(24c) 


Vl7^6 


- 7^6V2-l 


(24d) 


V27^6 


= 7^6VlV2 


(24e) 


ViCT 


= 0V2 


(24f) 




= a. 


(24g) 



Let us now determine the subgroup X,, of transforma- 
tions which are necessarily even. It evidently includes 
V^, V| and Tl\ (noted TZ-i). But there are more even 
transformations in this subgroup. Using Eq. 24c we find 
^V2^'Tia — ^'Re^Vi^V2J so =1. In the same way, using 
Eq. 24d, we get e^^ = 1- Thus X,, is generated by Vi, V2 
and TZ3. The algebraic relations in X^ are 



V1V2 = V2V1, 

7^3Vl = V27^3, 
7^3 = VlV27^3V2. 



(25a) 
(25b) 
(25c) 
(25d) 



As explained in Sec. IIIC, each of these relations gives 
a constraint on the gauge transformations associated to 
these generators. The Eqs.25 imply that for any site i: 

OvAV^'^)~OvAt) = PiTT, (26a) 
OnA^) + SnsiT^si) + OnAT^ji) = P2vr, (26b) 
Oils (i) - iV2'i) - Ov, (») = Ps^r, (26c) 

+0v.(V27^^^)-07^3(*) = P47r, (26d) 

where pi to p4 can take either the value or 1 (the equa- 
tions are written modulo 27r). We note [x] the integer 
part of X and x* — x — [x] (0 < x* < 1). By partially 
fixing the gauge, we can impose 

dvi{xi,yi) = 

Through a gauge transformation G of argument Oq, the 
6x of a lattice transformation X becomes: 



and the algebraic PSG is transformed in an other ele- 
ment of its equivalence class. Using the following gauge 
transformations: 



G4 



{x,y) 
{x,y) 



nx, 
Try, 



we see that a change of ps or p4 is a gauge transformation, 
so we can set them to zero. Solving the set of equations 26 
leads to: 



^Vi (x, y) 
&V2 i^: y) 





PlTTX 



{X,y) = PlTTX y - 



X + 1 

2 



(30a) 
(30b) 

gnd^*,y*),{iOc) 



with a supplementary constraint that can only be treated 
when the spin lattice is defined: 



9nA^*^y*) + ausii-y)*, (x-y)*) 
+gndiy - A-^Y) 



P27r. (31) 



This constraint only depends on the coordinates of the 
sites in a unit cell (x* and y*). 

Eqs. 30 and 31 define the chiral algebraic PSG on the 
triangular Bravais lattice. The full determination of the 
WS Antsatze requires precise definition of the spin lattice 
(triangular, honeycomb (m = 2) or kagome (m = 3)) and 
on the number of interactions included in the MF Hamil- 
tonian (first neighbor only or first and second neighbor; 
A and B parameters, or A only...). The case of the trian- 
gular lattice (to = 1) with nearest neighbor interactions 
and A and B MF parameters is described in the next 
subsection. 



V. STRICTLY AND WEAKLY SYMMETRIC 
ANSATZE ON THE TRIANGULAR LATTICE 
WITH FIRST NEIGHBOR INTERACTIONS 



Construction of WS Ansatze on the triangular 
lattice 



The triangular lattice has a single site per unit cell and 
the values of x* and y* are the coordinates of this site in 
a unit cell, say (0,0). Eq. 31 simplifies into: 



657^3(0,0) =0. 



(32) 



The solutions are g-jiaiOjO) = kn/S, with k integer. Be- 
cause the IGG is Z2, only the three values k = —1,0, 1 
lead to physically different Ansatze. 

Finally, we have 6 distinct algebraic PSG's for the re- 
duced set of symmetries Xe- They are characterised by 
two integers pi = 0, 1 and k — —1, 0, 1 and defined by: 



hi {x, y) 
?V2 (a;, y) 





PlTTX 



{i)-^eGii) + 9xii)-eGix-'i). 



(28) 



X + 1 



(33a) 
(33b) 

(33c) 
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FIG. 3: (Color online) Ansatze respecting the symmetries 
on the triangular lattice. All arrows carry Bij parameters 
of modulus B\ and of argument <j)Bi and Aij parameters of 
modulus A\ and of argument on red arrows (choice of the 
gauge), 2fc7r/3 on blue ones and 4fc7r/3 on green ones. On 
dashed arrows Aij and Bij take an extra pivr phase. 



Now, we have to find all the Ansatze compatible with 
these PSG's. The first useful insight is to count the 
number of independent bonds. Here, one can obtain any 
bond from any other by a series of rotations and trans- 
lations {i.e., elements of Xe)- Thus, if we fix the value 
of Aij and Bij on a bond ij, we can deduce all other 
bond parameters from the PSG. Note that Aij can be 
chosen real by using the gauge freedom. The value of all 
bond parameters are represented on Fig. 3 as a function 
of their value on the reference bond. The unit cell of 
the Ansatz contains up to two sites because pi may be 
non-zero. 

From now on we can forget about the PSG construc- 
tion and only retain the definition of the Ansatz given by 
Fig. 3 and its minimal set of parameters: two integers pi 
and fc, two modulus Ai and Bi, and one argument 4>Bi- 

Until now, we have only considered the subgroup X^. 
and we have looked for Ansatze strictly respecting these 
symmetries. We now want to consider all symmetries 
in X, but the symmetries in Xq will be obeyed mod- 
ulo an eventual time-reversal symmetry. This requires 
supplementary conditions on the Ansatze of Fig. 3. As 
explained in Sec. IV B, the transformations of Xo im- 
ply relations between the modulus and the arguments of 
the Ansatz. Since we are in a very simple case, where 
all bonds are equivalent in X^, no extra relation on the 
modulus can be extracted from Xq. However, some con- 
ditions can be found by examining how the the fluxes 
Arg(^ijy^*]^^ki-4[* ) on an elementary rhomboedron and 
Arg(^ijSjk.Aj^;) on an elementary triangle transform with 
TZq and a. Assuming that neither Ai nor Bi are zero we 
find: 



2A:^(l-eKj/3 = 

2fc7r(l + e^)/3 = 

(1 -I- eTZ^)(j>Bl = PlTT 

(1 - ecr)(f)Bi = PlTT 



(34a) 
(34b) 
(34c) 
(34d) 



For each set {e-ji^^ea), the compatible Ansatze are thus 
limited to: 

i) (e7^g,e^) = (1, 1): fc = 0, pi = and = or tt, 

ii) (e-R,,, Eo-) = (-1, -1): fc = 0, pi = and (/fs^ = or tt, 
in) {ens,f-c) = (1, -1): ipBi = Pi7r/2 or tt -|-pi7r/2, 

iv) i^-Rsi^o) = (^Ij !)• = 0, pi = and no constraint 
on cj)B^ ■ 

A couple (cTjej, Ecr) docs not characterize an Ansatz. A 
given Ansatz, can be found for several couples of parities. 
For example, the Ansatze obtained for (e^gje^) = (1, 1) 
are also present for all other {e-ji^^ea). Indeed as their 
MF parameters are real, they are not sensitive to time 
reversal and any t-ji^, can be chosen. From the clas- 
sical point of view, these Ansatze describe coplanar spin 
configurations, which are invariant under a global spin 
fiip followed by a TT rotation around an axis perpendicu- 
lar to the spin plan. 

Finally, there are nine different WS Ansatze families, 
given in Table I. We now conclude this section by a series 
of remarks concerning the solutions we have obtained: 

i) The number of WS Ansatze families is larger than 
the number of algebraic PSG of X^ , because the 
operators in Xq can act in different ways on the 
Ansatze. 

ii) Amongst these 9 Ansatze families, only the two first 
are non chiral, and the 6 others are TRSB Ansatze 
(by applying T, fc = 1 is changed to k — ~1 and 
(pB-^ to —ipBi)- The 6 families obtained by ap lying 
T are not listed here. 

iii) These solutions are called families as the moduli 
Ai and Bi can vary continuously without modi- 
fying the symmetries. The third Ansatz has no 
fixed value for <j)Bi and includes the first and sec- 
ond Ansatze families (they are kept as distinct as 
they are non chiral). 

iv) The fluxes of these Ansatze are easily calculated 
using Fig. 3. 

v) The detailed list of compatible Ansatze depends on 
the choice of the mean-field parameters (here, non 
zero Aij and Bij on first neighbor bonds) as we 
explain in App: C by contrasting these results to 
those of Wang et al. on the same lattice. 



B. Condensation of the WS Ansatze: the missing 
tetrahedral state 

The SBMFT has already been used to study the anti- 
ferromagnetic Heisenberg first-neighbor Hamiltonian on 
the triangular lattice with the Aij-O'abf decoupling^'-' 
(Eq. 5c) or with both Aij and S^^'' (Eq. 5a). The clas- 
sical limit of this model gives the well known three sub- 
lattice Neel order with coplanar spins at angles of 120°. 
The bond parameters obtained from this classical order 
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Ansatz number 


Pi 


k 


<t>Bi 


1 









2 







TV 


3 







any 


4 




1 





5 






TV 


6 







Tv/2 


7 


1 


3tv/2 


8 




1 


Tv/2 


9 




3tv/2 



TABLE I: The nine weakly symmetric Ansatze families on the 
triangular lattice, with the notations of Fig. 3. The modulus 
Ai and Bi axe not constrained although supposed non zero. 





F 


coplanar 


tetra 


Pi 








1 


k 







1 


Sr 


? 


? 


1 


Scr 


? 


? 


-1 







Va* 

4 


1 • 

76 


Bi 


1 * 

2 


1 • 
4 


1 * 
yi2 


4>Bi 





TV 


7r 

"2 



TABLE II: Values of the parameters of Fig. 3 for Ansatz fami- 
lies related to regular classical states on the triangular lattice. 
The states are designed by F for ferromagnetic, coplanar is 
the \/3 X \/3 state and tetra for tetrahedral. These states are 
described in more details in . The interrogation points mean 
that the two values e = ±1 are possible (coplanar or colinear 
state). The * means that the parameter value is free, we give 
its value in the fully magnetized state (jm| = 1). 



(see Eq. 21) lead to a strictly symmetric Ansatz (no need 
to break T: we can chose to fix (e-T^gje^) = (1,1)) with 
Pi = 0, fc = and — ti"- We note that all MF param- 
eters are real in this gauge choice (is it always possible to 
do so for coplanar states). In this case the restriction to 
real bond parameters did not prevent to obtain the true 
MF ground state. 

The tetrahedral state (Fig. 1) is the unique GS of 
the multi-spin exchange Hamiltonian in a large range of 
parameters:''''^" 



in) 



J A 



(ijkl) 



{P{ijkl) + P{ilkj)), (35) 



where the second sum runs on every elementary rhom- 
boedra and P(ijki) is a cyclic permutation of the spins 
and J4 > and | < ^ < 1. Moreover, it is one of the 
GS's of a Heisenberg Hamiltonian with first and second 



neighbor interactions 

iij) ((ij)) 



(36) 



for i < a < 1. In the later situation the GS is how- 

o 

ever degenerate and fluctuations (order by disorder) fa- 
vors coUinear orders. 

The bond parameters obtained from this classical 
order (Eq. 21) lead to the weakly symmetric Ansatz 
((eKe^Ccr) = (1,-1)) with = 1, fc = 1 and (/)Bi = 7r/2 
(or opposite k and (/>Bj for the opposite chirality). The 
previous SBMFT studies of the ring exchange model 
(Eq. 35) have been limited to real parameters and it 
would be interesting to perform a systematic search for 
a possible chiral MF ground-state. If the chiral Ansatz 
indeed turns out to have the lowest energy - as suggested 
by its classical limit - then the spin-^ might be a chiral 
SL since exact diagonalisations"'''''^" have shown the ab- 
sence of Neel long range order in some parameter range. 



VI. FLUXES 

We have already given a brief definition of the fiuxes 
in Sec. Ill A; in this section we will enlarge this definition 
and comment on the physical meaning of the various loop 
operators (local and non local) that can be defined on a 
lattice. 

The gauge invariance of a product of A^j, Aj^-, Bij 

and operators on a closed contour requires two con- 
ditions: (i) each site i appears in an even number of 
terms, (ii) The set of operators containing a site i can 
be organized into pairs such as the product of each 
pair is invariant by a local gauge transformation on site 
i (for example Aji and Bik). Such a gauge- invariant 
operator is the analog of a Wilson loop operator in 
gauge theory and the complex argument of its expec- 
tation value is called a flux. Aig{AijA^-f,...AimA\^^)^ 



Ajcg{BijBjk ■ ■ ■ Bii) are examples of fluxes with only Aij 
or Bij operators, but it is possible to mix both as for 
example in Arg(A„A]^i?^;l,,„A|„,). In SBMFT we ap- 
proximate these averages of products by the product of 
the averages (this can be formally justified in the N 00 
limit). For example: {B^^Bjk ■ ■ ■ Bu) BijBjk ■ ■ ■ Bu. 

There is an infinite set of non-independent fluxes. "^'^ A 
method to determine the number of independent fluxes 
for a given set of non zero Aij and Bij is given in App. E. 
To characterize a given Ansatz, we can limit ourselves to 
the minimal set of independent parameters that define 
unequivocally its equivalence class: essentially the non- 
zero bond field modulus and a minimal set of fluxes. 

The first insight on the physical meaning of the fluxes 
is given in the classical limit (Sec. VIA), where they are 
simple geometric quantities related to the orientation of 
the spins. Then, we come back to the quantum case and 
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express the fluxes, which are physical quantities, with the 
exclusive use of spin operators (Sec. VI B). 

A. Definition and physical meaning in the classical 

limit 

We first concentrate on the mean-field flux formed by 
products of Bij parameters. In the classical limit, the fiux 
of Bij around a loop ijk . . .1: Arg(Sij;Bjk . ■ . B\{) is related 
to the solid angle associated to the contour described by 
the spins on the Bloch sphere. We give here a simplified 
formulation of the calculation given in Ref. Ki. Let us 
suppose that the direction of the magnetization (with a 
modulus fixed to 1) evolves slowly along the loop and use 
the gauge of Eq. 20, but in spherical coordinates: 




Arg(Si*) ~ S(l - cos^O^^y^. (38) 

This last quantity (to first order in the variation of the 
spin) is the half of the solid angle between the three di- 
rections defined by the z axis and the spins at site i and 
J. By summing such quantities around a closed contour, 
we obtain the half of the solid angle spanned by the spins 
along the loop. This illustrates the gauge dependence of 
a single S*^-: by a gauge transformation we change the 
direction of the z axis and thus Arg(,B^), but the total 
solid angle of the closed loop is independent of the choice 
of the z. 

In a similar approach the flux 
Arg(Aj(-^*j^) . . . Ain(--4*ji) is associated to the 
half of the solid angle defined by the spins along the 
loop, but after flipping one spin every two sites (the j 
spin for Aij, the i for —A*j). The — I's present in the 
above expression have their importance as they can lead 
to a final difference of tt. 

For more complicated fluxes mixing Aij and Bij pa- 
rameters, we flip one spin every two sites on Aij and 
A^j bonds (as previously), we flip all of them for Bij, 
and none for B*j. The flux is then half the solid angle 
associated to these modified spin directions. 

We can now reformulate the previously discussed re- 
lation between chirality and fluxes. If a classical state 
is chiral, it has non trivial fluxes on contours where the 
spins are non coplanar. If the corresponding MF param- 
eters are non zero, we then have found a loop with a 
non-trivial flux and whatever the gauge choice, at least 
one MF parameter has to be complex. Now, if a state is 
coplanar, then all fluxes are trivial and in a gauge where 
the spin plane is xz, all MF parameters are real. 

In the tetrahedral state described on Fig. 1, the fiux of 
the Aij around a small rhomboedron is ±7r/3 and the flux 
of the Bij around a small triangle is ±7r/2 (depending on 
the choice k = ±1, see Sec. VB). 



B. Fluxes in quantum models 

In the quantum realm, the fluxes can no longer be ex- 
pressed in term of solid angles. But as we have already 
noted, Wilson loop operators are gauge invariant quan- 
tities and as such, they are physical observables and can 
be expressed in terms of the spin operators. 

1. Spin- 1/2 formulas 

To simplify we will start by imposing that the con- 
straint is strictly verified for S = |, so there is exactly 
one boson per site. We have noted that in the classical 
limit, the scalar chiralities are associated to the fluxes. 
In the quantum case, we can express the flux operators 
in term of permutation operators, generalizing some re- 
sults of Ref. 7. The operator that transports the spins at 
sites 1, 2, 3 to sites 2, 3, 1 is the permutation noted P(i23). 
We recall that the permutation operator of spins between 
two sites can be written as: 

P(,,) = i -f 2S, • S, (39) 

This straightforwardly implies that the fiux of the Bij 
operators is 

: BlBl...Bi, := ^P(i2..„) (40) 

The formula for the fiux of the Aij operators is more 
involved. It reads 

: yl^2^23^34 ■ ■ ■ ^2nl ^P{l2..2n){^ - P{23)) 

(l-P(45))---(l-i5(2„l)). 

To prove this last assertion, we first note that ^'''^ 
is the projector on the singlet state of the two spins i 
and j. We then verify this equality in the basis of states 
^7=i'^2i,2i+i, where ipij are eigen vectors of P(ij)- In 
the case where at least one bond is in a symmetric state 
(triplet), both sides of Eq. 41 are zero. The final step 
is simply to check that the relation holds for the state 
which is a product of singlets. 

2. Fluxes in quantum spin S models 

For S > 1/2, Eq. 39 is no more valid and Eqs. 40 and 
41 are not more valid either. But we can still replace the 
on-site number of bosons by 25* and obtain an expression 
depending only on the spin operators. The expression of 
the product of four Aij operators is: 

8 : ll2l23^L^41 := (Si.S2)(S3.S4) + (S2.S3)(S4.Si) 
-(Si.S3)(S2.S4) + 52(81.83 + S2.S4 - S1.82 

-82.S3 - S3.S4 - S4.S1) + 54 + z5(84.(Si X 82) 

-Si.(S2 X 83) + 82.(83 X 84) - 83.(84 X 81)) 



13 



The expression of the product of three Bij operators is: 

4 : bI^bIMi ■■= 5(Si.S2 + S2.S3 + S3.S1) 

+8^-181.(82x83) ^ ' 

3. Fluxes m SBMFT 

In a state where on-site number of bosons is not strictly 
conserved, the previous expressions become a bit more 
comphcated. The number operators can no longer be 
replaced by 25', and we have for example: 

4 : bI^bI^bI := ingSi.Sa + iriiSa-Sg + ^^283.81 

+ ^lf^-iSi.(S2 X S3). 

(43) 

C. Finite size calculations lattice symmetries and 
non local fluxes 

For simple lattices as the square or triangular lattice, 
we can solve analytically the MF Hamiltonian -ffMF of 
Eq. 9 directly in the thermodynamical limit. But in most 
cases, we have to solve numerically the self-consistency 
conditions on finite lattices. 

To use the chiral PSG's on a finite periodic lattice, 
we have to be cautious about symmetries. Indeed, all 
precautions have been taken so that the Ansatz (strictly 
or weakly) respects the lattice symmetries on an infinite 
lattice. But we have to verify that the finite periodic 
lattice has the same symmetry group as the infinite one. 
This verification is quite usual for local properties, but 
is more subtle for non-local ones and can be most easily 
understood in term of fluxes on large non-local loops. 

PSG's impose that fluxes on local loops are preserved 
by lattice symmetries (or sent to their opposite in the 
case of a chiral state) . But some additional care has been 
taken concerning loops which are topologically non triv- 
ial (cannot be shrunk to a point by a succession of local 
deformations). These loops which "winds" through the 
boundary conditions do not exist on the infinite lattice. 
For a symmetric Ansatz to remain symmetric on a finite 
periodic lattice, we have to verify that the fluxes associ- 
ated to these topologically non trivial loops also respect 
the lattice symmetries. The way to treat the problem of 
the non-local loops is detailed in App. F, together with 
several ways of understanding their meaning. 

VII. CONCLUSION 

In this paper we have extended the PSG construction 
to include time-reversal-symmetry-breaking states with 
the SBMFT. These TRSB phases that we describe gener- 
ically as chiral , can also break one or many discrete sym- 
metries of the lattice (in the triangular example either a 
or 72-6 ). Using this constructive method we have built 



all the SS and WS Ansatze with two MF parameters on 
the triangular lattice. All the regular 0(3) magnetically 
ordered phases can be obtained from these Ansatze by 
spinon-condensation (the others have no regular classical 
limit). The TRSB Ansatze have, when they condense, 
non-planar magnetic order and non-zero scalar chirali- 
ties. 

The TRSB SL have short range spin-spin correlations 
but non trivial fluxes on various loops. The simplest of 
these fluxes are related to the imaginary part of the per- 
mutation operator of three spins, that is directly related 
to their scalar chirality. In some cases the time-reversal 
symmetry breaking fluxes might be more complex, as ex- 
plained in section VI and illustrated in Appendix D 2 for 
the kagome lattice. These various fluxes have been ini- 
tially defined within the MF Schwinger boson approach 
but VI has shown how these gauge invariant quantities 
can be expressed in terms of spin operators, indepen- 
dently of any MF approximation. It should be noticed 
that in a TRSB SL fluxes other than those deduced from 
the Ansatz may be non zero and easier to compute. It is 
the case for example in the cubocl SL recently proposed 
for the nearest- neighbor Heisenberg model on the kagome 
lattice. '^ The flux of the A bond operators around the 
hexagons can be expressed in terms of spin permutation 
operators but it is relatively involved (Eq. 41) and has 
not yet been computed numerically. In fact, in that phase 
(at least at the MF level), there are simpler fluxes which 
are non zero, as for example the triple product of 2nd 
neighbor spins around hexagons, or the triple product of 
three consecutive spins on an hexagon. 

In Sec. VB, another TRSB Ansatz was discussed in 
relation to the ring exchange model on the triangular 
lattice. 

In spite of short range spin-spin correlations the TRSB 
SL have some local order parameter associated to the 
fluxes. The finite temperature broken symmetries being 
discrete symmetries, there are no Goldstone modes and 
these chiral phases should survive thermal fluctuations 
in 2D. The phase transition associated to the restoration 
of the chiral symmetry has been studied in some classi- 
cal spin models. ^^''^"''^^ In spite of the Ising like character 
of the order parameter, the phase transition was shown 
to be weekly first order due to interplay of vortices in 
the magnetic texture with domain walls of the chirality. 
It has been shown within the SBMFT framework in the 
cubocl phase that thermal fluctuations tend to do ex- 
pel the chiral fluxes (favour coplanar correlations) but 
a more complete study (beyond MF) of the flnite tem- 
perature properties of a TRSB SL would be required to 
understand the specific properties of the chiral transition 
in these systems. 

Finally, it would be useful to clarify the "topological" 
differences (entanglement, degeneracy, edge modes, ...) 
between the present chiral SL described in the SBMFT 
framework with the chiral SL wave-functions related to 
fractional quantum Hall states (such as the Kalmeyer- 
Laughlin state or that of Yang, Warman and Girvin for 
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instances) , as well as the difference with conventional (T- 
symmetric) Z2 liquids. It would also be very interesting 
to analyze qualitatively the effects of (gauge) fluctuations 
in the present chiral SL. 



Appendix A: The Bogoliubov transformation 

This appendix explains how to obtain the eigenmodes 
of Eq. 9. New bosonic operators, components of (jj, are 
created by linear combinations of the components of (f) to 
obtain a new diagonal matrix M. For the Hamiltonian to 
possess a GS (spectrum bounded from below), the diag- 
onal elements (wi, . . .^2^, ) must all be positive or null. 
This transformation is called the Bogoliubov transforma- 
tion and is generally well documented (see for instance 
Ref. IG) when the size of the matrix M is 2 x 2 (the 
transformation can then be done analytically), but more 
rarely for larger sizes (where numerical calculations are 
sometimes required) . When periodic Ansatze are consid- 
ered, a Fourier transform can block-diagonalize M, with 
blocks of size 2m x 2m, with m the number of sites in the 
unit cell. As soon as m > 1, Bogoliubov transformation 
of matrices larger than 2x2 are needed. 

Note that the choice of an Ansatz whithout any Aij 
parameters (for example using Eq. 5b) simplifies con- 
siderably the Bogoliubov transformation since the total 
number of boson is conserved and M is block diagonal 
with two blocks of size N^. The transformation reduces 
to the diagonalization of each block by a unitary matrix. 
The new bosons bi^ are then linear combinations of the 
old bia, without any component. The vacuum of the 
new bosons is the same vacuum as for the old bosons. To 
respect the constraint on the boson number, we have to 
create a Bose condensate (see Sec. IV A), which implies 
long-range magnetic order. This proves that the Aij pa- 
rameters are necessary to obtain SL. 

Here we describe the general method for the cases 
where M can have an arbitrary size, as explained in de- 
tails in by Colpa. The 2Ns x 2Ns matrix P defined 
such that (j) = Pcj) is called the transformation matrix. 
Let us look at the conditions P should satisfy. The most 
evident is that M must be diagonal, which gives a first 
constraint. The second one is that the Ng first compo- 
nents of (j) must be annihilation operators and the last 
Ns, creation operators. This gives a constraint on their 
commutation relations. The two resulting conditions are 



P^JP ^ J 



(Al) 



where J is the 2Ns x 2Ns diagonal matrix with coeffi- 
cients — 1 for the Ns first terms and 1 for the last Ns 
elements {Jij = [(/)J,0j]). The second constraint makes 
the Bogoliubov transformation different from a diagonal- 
ization (where J would be the identity matrix). It is 
sometimes called a para-diagonalization. 

Here we just recall the main steps of the algorithm ''^ 
to solve these equations: 



• Verify that M is definite positive. It ensures that 
the GS is unique (in some cases where M has zero 
eigenvalues, the GS exists but is not the unique). 

• Find a complex upper-triangular square matrix K 
such as Al = K^K (Cholesky decomposition of M). 

• Find a unitary matrix U such as i = W KJK^U 
is diagonal with it first Ns coefficients positive and 
the other negative (usual diagonalization of a her- 
mitian matrix). 

• The solution is M = JL and P = R-^UM^I^. 

Using the rotational invariance, we deduce that the Ns 
first coefficients (wi, . . . of M are the same as the 

Ns last (maybe differently ordered). The energy of the 
MF Hamiltonian GS writes: 



1=1 



(A2) 



and its elementary excitations are free bosonic spinous 
with energies (wi, . . . , i^jv^) and spin 1/2, from which we 
can get the free energy at any temperature. We are now 
able to look for solutions of Eq. 11 and 12, i.e. the sta- 
tionary points of the free energy with respect to the MF 
parameters and with respect to the Lagrange multipliers. 



Appendix B: Bounds on self-consistent values of the 
MF parameters in SBMFT 

The moduli \Aij\ and \Bij\ are a priori unconstrained 



real numbers in SBMFT. We prove here that in a self- 
consistent Ansatze, their moduli cannot exceed an upped 
bound: \Aij\ < ^^y^ and \Bij\ < ^. These inequalities 
considerably restrict the domain to explore and facilitate 
the numerical search for solutions. 

Let be any normahsed bosonic state. We denote 
by (O) the the expectation value of an operator in this 
state. Whatever the operators u and v we have 



l(™>l < 



Applying it to Aij and Bij, we obtain 



(Bl) 



1(4)1 < ^^'+7+^\ 1(4)1 <^^. (B2) 

We now take {(p) as the GS of Hmf (Eq- 9) for some 
Ansatz. If the chemical potential is adjusted, (n^) = k 
on every lattice site. In the case of self consistent param- 
eters, \Aij\ = 1(4)1' = 1(4)1 and Eq. B2 leads 
to 



(B3) 
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Appendix C: The strange classical limit of the tt flux 
Ansatz of Wang and Vishwanath ' 

Wang and Vishwanath"' explored all the strictly sym- 
metric Ansatze {(e-izg,e^) — (1,1)) with the Aij decou- 
pling for first neighbor Heisenberg interactions. They 
found two Ansatze. The first one is characterized by 
a flux Aig{Ai}A^^^A]iiAl^) = around a rhomboedra for 
(pi = 0), giving the 3 sublattice Neel order in the classical 
limit. The second one has a flux tt rhomboedra {pi — 1). 
The Aij parameters they used are those obtained from 
Fig. 3 with the corresponding value of pi and fc = 0. 

Comparing this to our result for the SS Ansatze, we 
may wonder why do they obtained two possibilities for 
Pi (0 or 1) whereas we found that pi = was the only 
solution for (e^g, e^) = (1, !)■ The difference comes from 
the the absence of Bij parameter in their MF approach. 
The complex phase of Bij is then ill defined and only 
the first of our two constraints (Eq. 34) remains. They 
thus impose k = 0, but nothing on pi. In fact, as this 
situation is the limit i?i of none of the SS cases we 
have explored in Sec. V, it appears that the 7r-flux Ansatz 
is unstable with respect to the introduction of Bij. In 
other words, any non-zero value of Bij will break at least 
one lattice symmetry. The two WS Ansatze described in 
the 6th and 7th lines of Tab. I corresdpond to this limit. 

The nature of the spinon condensation in the 7r-flux 
Ansatz could not be completely clarified in Ref. ■. Our 
understanding is that it is not consistent to impose Bij — 
to describe ordered states on a frustrated lattice. The 
only way to have \Bij\'^ = classically is indeed to have 
anti parallel spins on all bonds, which is not possible on 
the triangular lattice. 



Appendix D: Weakly symmetric Ansatze on some 
usual lattices 

1. Lattices with a square Bravais lattices 

The first step is to find all chiral algebraic PSG's. We 
choose the most general case IGG^ Z2 and we suppose 
that Hq respects all the lattice symmetries whose gen- 
erators are described in Fig. 2 (right). The coordinates 
(x, y) of a point are given in the basis of the the trans- 
lation vectors Vi, V2 and the action of the generators on 
the coordinates are 

Vi : {x,y) ^ ix + l,y), (Dla) 

V2:ix,y) ^ ix,y + l), (Dlb) 

TZi : {x,y) {~y,x), (Die) 

cr : {x,y) (y,x). (Did) 



The algebraic relations between them are 



V1V2 


= V2V1 


(D2a) 


V27^4 


= n^Vi 


(D2b) 




= I 


(D2c) 


Vl7^4V2 


= 7^4 


(D2d) 


Vicr 


= aV2 


(D2e) 




— a 


(D2f) 


a' 


= I 


(D2g) 



To our knowledge, even the non-chiral algebraic PSG's 
have not been derived previously. Here, we directly de- 
rive the chiral ones. From Eq. D2, we deduce that the re- 
duced set of symmetries is generated by , V| and TZi 
(noted V'l, V'2 and 1^2) ■ Moreover, we find that evi — e-Vi- 
An Ansatz is characterised by the parities (ey^ , £7?,, Co-)- 

The algebraic relations between these generators are 



V[V^ - V^V(, (D3a) 

nl = /, (D3b) 

V{R2V[ = 7^2, (D3c) 

V^7^2V^ = 7^2. (D3d) 



As explained in Sec. IIIC, each of these relations gives 
a constraint on the gauge transformations associated to 
the generators. The constraint from Eq. D3 are then, for 
all i: 

^v^(Vr'i)-0v^(i) = Pi7r(D4a) 
OnS) + SnAT^2i) = P27r(D4b) 
en,{V[-\) - OnA^ = P37r(D4c) 
Ov, (V^z) + {n2i) + On, (*) - On, (^) = P47r(D4d) 

where pi , • • • , P4 can take either the value or 1 (the 
equations are written modulo 27r). We note [x\ the in- 
teger part of a;/2 and x* — x — 2[x\ (0 < a;* < 2). By 
partially fixing the gauge, we impose 

dv[{xi,yi) = 0, 
Qv:^{x*i,yi) = piTTx*. 

Contrarily to the triangular lattice, no gauge transforma- 
tion can here be used to get rid of some pi . 
Solving the previous equations D4 leads us to 

evAx,y) = 0, (D6a) 

Sv2{x,y) = PiTTX, (D6b) 

^n2{x,y) = P3TTX + p4TTy + gn2{x* ,y*)- (D6c) 

with a complicated supplementary constraint that can be 
treated only when the lattice is more precisely defined: 

g^Jx*, y*) + <?K. ihxy, i-y)*) = P27T. (D7) 

This constraint only depends on the coordinates of the 
sites in a 2 X 2 unit cell (x* and y*), so it gives at most 
4m independent constraints. 

These general algebraic PSG's can then be used to 
find the weakly symmetric Ansatze on any lattice with 
a square Bravais lattice (for example: the square, the 
Shastry-Sutherland lattice,. . . ). 



16 



Weakly symmetric Ansatze on the kagome 
lattice 



The Bravais lattice of the kagome lattice is triangular, 
so we use the algebreaic PSG's determined in Sec. IV C. 
The unit cell contains three sites. We choose to place the 
origin of the frame at the center of an hexagon and the 
coordinates of the sites in a unit cell are (|, 0), (0, and 
f- -) 

Since the sites have non-integer coordinates, it is con- 
venient to transform the Eqs. 30 using the following 
gauge transformation (see Eq. 28): 



Gi : {x,y) - 
The new algebraic PSG is 


Pi7r[x] 

+9n-Ax* ,y*), 



■piiryx 



On^x^y) 



(D8) 



(D9a) 
(D9b) 



N + 1 



[y* 



(D9c) 



Even if it seems more complicate than Eq. 30, it avoids 
some pi7r/2 and simplify the future Ansatze. This gauge 
transformation is equivalent to a different initial choice 
of evM.y^) in Eq. 27: 



(DIO) 



Under the effect of G, Eq. 31 is modified and gives the 
constraint 



5K3 



1 



I + 9n, ( 0, ^ I + gn. 



1 1 

2' 2 



= (P2+Pi)7r. 



Using the following gauge transformations. 



G2 

G3 
G4 



{x^y) — ax*, 
{x,y) -> by*, 
{x,y) [y* - x*]Tr, 



(Dlla) 
(Dllb) 
(Dllc) 



with a and b real numbers, we can set 177^3 = 0. Finally, 
we have two distinct algebraic PSG's for the reduced set 
of symmetries. They are characterised by pi = ±1 and 
defined by Eq. D9 with g-j^^ = 0. 

Now, we have to find all Ansatze compatible with these 
PSG's. We limit ourselves to first neighbor parameters, 
but this procedure is easily generalized to further neigh- 
bors. Two bonds are needed to generate the whole lattice 
by action the reduced symmetries: one blue bond and 
one of red bond of Fig. 4. The values of Aij and Bij on 
these two reference bonds allows to obtain all the other 
bond parameters by using the PSG. Note that Aij can 
be chosen real for say the reference blue bond by using 
the gauge freedom. The values of all bond parameters 
are represented on Fig. 4 as a function of their value on 
the reference bond. The general unit cell of the parame- 
ters contains six sites because of the possibly non-zero pi . 




FIG. 4: (Color online) Ansatze respecting the Xe symmetries 
on the kagome lattice. Blue arrows carry Bij parameters of 
modulus Bi and of argument (jjBi and Aij parameters of mod- 
ulus Ai and of argument 0. Red arrows are for modulus B[ 
and A'l and arguments (I>b[ Bjcid . On dashed arrows Aij 
and Bij take an extra Pitt phase. 



Note that the simplicity of the Ansatz of Fig. 4, where 
phases differ only by tt between two bonds of the same 
color is a consequence of the choice of Eq. DIO. 

Finally, we can forget all about the PSG and only re- 
tain the parameters needed to completely describe an 
Ansatz together with Fig. 4. These parameters consist 
in the integer pi , the modulus Ai , A[ , Bi and B[ , and 
the arguments c/)^'^, (pBi and (/j^j. 

We now want to consider all symmetries in Xq. As 
blue and red bonds are related through TZq, this implies 
relations between the modulus: Ai — A[ and Bi = B[. 
Relations on the phases are found by looking at the effect 
of Tie and a on the flux Aig{AijA*i^AkiAl^AmnA*^i) of an 
elementary bow tie and Arg(^ijSjk^j^i) of an elementary 
triangle (we suppose that neither Ai nor Bi are zero): 



(cKe + ea)0Si = 

(1 + = 



(D12a) 
(D12b) 
(D12c) 
(D12d) 



For each couple {efig,ea-), the number of compatible 
Ansatze is thus reduced: 

i) (eKe'^tr) = (1, 1): (f>Bi = or tt and (f)A[ = or 



= (1,-1): ^B, = 



and 



~ or TT, 



iii) {<^i16,<^'t) = (-1, 1): '/'Si = 

iv) (e-Re-e^r) = (-1,-1): (f>Bi = 



and = or TT, 
and (j)Bi = or TT. 



Finally, there are 20 different WS Ansatze families, 
given in Table III. Each regular states of Ref. 9 belongs to 
one of them: the 2nd for the q , the 6th for the \/3 x 
-v/S , the 17th for the octahedral, the 20th for the cubocl 
and the 14th for the cuboc2 state. The parameters of the 
fully magnetized states are calculated using Eq. 20. They 
are also used to determine the parities and pi. In the 
case where both parities are possible (coplanar states) , we 
chose the largest symmetries (e = 1) to fix the maximum 



17 



Ansatz number 


Vi 




<f>Bi 




1 











2 









TT 


3 








any 


4 








any 


5 












6 




TT 




TT 


7 








any 


8 








any 


9 




any 







10 






TT 


11 












1 


same 


as for pi = 


20 











TABLE III: The twenty weakly symmetric Ansatze families on 
the kagome lattice, with the notations of Fig. 4. The modulus 
A\ = A'l and Bi — B'l axe not constrained, except that they 
do not vanish. 





F 


g = 




oct 


cubocl 


cuboc2 


Pi 











1 


1 


1 


£r 


7 


7 


7 


1 


-1 


-1 




7 


7 


7 


-1 


-1 


1 







2 


^* 

2 


1 • 
%/2 


2 


1 • 

2 


<i>A\ 







TT 


TT 


TT — atan\/8 





Bi 


1* 


1 * 

2 


1 • 
2 


1 • 


1 • 

2 


2 







TT 


TT 


4 


TT 


atanA/2 — tt* 


<f>B' 





TT 


TT 


■/"Si 


TT 


-0Si 



TABLE IV: Values of the parameters of Fig. 4 for Ansatz fam- 
ilies related to regular classical states on the kagome lattice 
(the states are designed by F for ferromagnetic, oct for octa- 
hedral and cuboc for cuboctaedron order parameters. These 
states are described in more details in ). The modulus verify 
A\ = A'l and Bi — B[. The interrogation points mean that 
the two values e = ±1 are possible (coplanar state). The * 
mean that the parameter value is free, we give its value in the 
fully magnetized state (|in| = 1). 



number of parameters. They are indicated in Table IV. 
For the parameters which are not fixed, we give the values 
take in the classical limit. The self consistent parameters 
for finite S are different, but generally not far from the 
classical values. Thus, they can be used as a starting 
point in numerical optimizations. 

These calculations are easily generalised to further 
neighbors and have already been used for two studies 
on the kagome lattice^^'^ 



Appendix E: Number of independent fluxes on a 
lattice 



We suppose that we have a MP Hamiltonian with 
TT'A + nis non-zero bond parameters {{Aij}, {Bij}) {Aij 
and Aji count as only one parameter, and the same for 
Bij and Bji). As they are complex numbers, we need 
2n_4 4- 2ne self-consistent conditions to solve this MP 
problem. We already know that the solution is not unique 
and that two Ansatze related by a gauge transformation 
are equivalent. Thus, by fixing the gauge, we can de- 
crease the number of equations for the complexe phases. 
In fact, the number / of necessary arguments corresponds 
to the number of independent fluxes on the lattice. In 
this appendix, we describe a simple method to compute 
/ on a finite cluster. 

We define a rectangular matrix M of size {nj^^+nig)xNs 
{Ns is the number of sites), where each line characterizes 
a MP parameter, and is therefore associated to a pair 
of sites (ij). As for the column, they correspond to the 
lattice sites. The coefficients of a line are all zero except 
for the two entries at columns i and j. Both entries are 
equal 1 for an Aji-hond, whereas these entries are —1 and 
1 for a Bji-hond (which site is ± 1 has no importance). 
Then the result is: 



/ = n_4 + nis — rank(M). 



(El) 



The effect of a gauge transformation on the bond 
phases is obtained by multiplying M by the vector 
(^1, • • -^jVs)*- By definition, a product of bond param- 
eters defines a fiux if the sum of their complex phases 
is unchanged by a gauge transformation. It means that 
the sum of associated matrix lines is 0. As the complex 
conjugate of a bond parameter can be used, the weight of 
each line in the sum can be ±1. As we can imagine using 
several times the same parameter, the weight of each line 
in the sum can finally be any relative integer. So, the 
existence of a flux relating a set of parameters is equiva- 
lent to the existence of a vanishing linear combination of 
their lines. 

We can now give the proof of Eq. El by induction. The 
relation Eq. El is true for one parameter: 



M 




M 




(E2) 



We suppose now Eq. El true for + parameter and 
we add a parameter on a bond (possibly with a new site) . 

• If a new site is added, the matrix gains a column 
and a line with a 1 at their intersection, so the rank 
of M increases by 1 and / remains the same. As 
we can chose the gauge on the new site, the new 
parameter can be chosen real, and Eq. El remains 
true. 

• If no new site is added and there is new flux us- 
ing the new parameter, the new line is a linear 
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combination of previous lines, thus, rank(M) is un- 
changed and / increases by 1. 

• If no new site is added and no new flux exists using 
the new parameter, the new line can not be writ- 
ten as a linear combination of previous lines and / 
remains the same. 



Appendix F: Example of non local fluxes breaking 
the lattice symmetries 




We illustrate the possibility for an Ansatz to be in- 
compatible with a periodic lattice. The example we 
give is an Ansatz on a 12-site periodic triangular lattice 
(Fig. 5) that strictly respects the infinite lattice symme- 
tries (Fig. 2). Let us choose for simplicity the Ansatz 
with only first neighbor Aij MF parameters defined by 
fc = and pi — 1 (already discussed in the classical limit 
in App. C). This Ansatz is the simplest illustration of 
these symmetry issues, but they can be encountered for 
any other Ansatz (as will become clearer later on). 

Periodic boundary conditions defining a finite lattice 
are defined by the two vectors: Ci and £2- Two sites 
separated by an integer linear combination of the £1 and 
C2 vectors (Fig. 5) are identified as the same sites. The 
three loops £1, £2 and £3 are mapped onto each other by 
rotations and should therefore have the same fluxes in a 
WS or SS Ansatz. But here their values are for £2 and 
£3 and TT for ii . This is due to the fact that the unit cell of 
the Ansatz is twice the unit cell of the triangular lattice, 
and it introduces a distinction between the directions Ci 
and £2 ■ The TZ^ symmetry cannot be restored simply us- 
ing a gauge transformation. We see in the figure that the 
combination of the three loops (^"1 -1-^2 +^3) is a local loop, 
with trivial winding numbers. Thus the flux of this loop 
is fixed by the Ansatz and it is tt. If we do not change 
the local physical properties of the Ansatz, which we did 
not want to, the sum of the three fluxes should remain 
equal to TT (modulo 27r), and the only way out is to have 
a TT flux on the three non local loops. This can be done 
by choosing a specific non local contour (here the green 
£1 contour for example) and adding an extra phase tt to 
all MF bond parameters crossing this contour. Trans- 
formation does not affect the local (contractible loops) 
fluxes since they always contain an even number of al- 
tered parameters. But the fluxes associated to £2 and £3 
acquire an extra phase factor. The three fluxes around 
£1 , £2 and £3 are all equal to tt and the symmetries of the 
infinite lattice are now respected for this finite periodic 
lattice. 

We can see this modification as a change in the bound- 
ary conditions (BC's), from periodic in both directions 

to periodic in the Ci direction and antiperiodic in the C2 
direction: 



FIG. 5; (Color online) Triangular lattice with 12 sites (filled 
circles) and periodicity Ci and £.2, respecting all the lattice 
symmetries of Fig. 2. Three non local loops £1, £2 and £3 
are drawn in dashed arrows. The green dotted lines cross the 
bonds where the MF parameters are multiplied by —1. 



where ni and n2 are arbitrary integers. This changes the 
set of allowed wave vectors k from 

k • £2 = 

to 

k-£i = 

k • £2 = TT 

The wave vectors of the 12-site lattice before and after 
the transformation are drawn in Fig. 6. The spinon dis- 
persion computed in the thermodynamic limit has two 
minima (dark red). Periodic boundary conditions for this 
12-site sample present evident drawbacks: the pattern of 
allowed wave vectors (blue points in Fig. 6) does not re- 
spect the 7^.3 symmetry of the spinon dispersion and the 
minimum of the spinon dispersion is not reached in the 
12 sites samples with this PBC. We could hastily have 
supposed that single-spinon states are not physical ex- 
citations and as such they do not have to respect the 
lattice symmetries. But this statement is incorrect. The 
vacuum of spinous calculated from the set of wave vec- 
tors obtained from periodic boundary conditions is itself 
distorted and so are any physical quantities as for exam- 
ple spin-spin correlations that are calculated from this 
input. On the contrary the modified BC's restore the 
TZq symmetry of the pattern of authorized wave vectors 
around the spinon minima. 

This can also be understood in a different way. The 
periodic or antiperiodic boundary conditions define the 
4 topological sectors on the torus. To go from one sec- 
tor to an other, we create two visons, '' ' move one of them 
around the lattice and annihilate them again. It is equiv- 
alent to the sign change of the MF parameters along this 
loop. The present discussion shows that for the 12-sites 
sample PBC do not define the (0,0) topological sector of 
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the model and we have to go to the APBC to describe 
this fully symmetric sector. In classical terms a change 
of TT in a flux around a loop corresponds to a rotation 
of 2tt of the spin orientations, thus to a Z2 vortex. ''''' 
Translated to the classical limit, the previously used pe- 
riodic BC's thus correspond to a twist of 2tt around the 
lattice, that is at the existence of a non trivial vorticity. 



FIG. 6: (Color online) The hexagon is the Brillouin zone of the 
triangular lattice, the rectangle, this of the Ansatz. The wave 
vectors of the 12 site lattice with periodic boundary conditions 
(BC's) are the blue ones, whereas they are the green ones for 
periodic EC in the £1 direction and antiperiodic in the £2 
direction. The background intensity is the value of the spinon 
energy (dark for minima). 
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tuations, leading eventually to a non chiral SL phase. But 
this quantum phase transition has no reason to be concomi- 
tant with the opening of a (spinon) gap and the appearance 
of a SL, as can be seen in Ref. 12. 

SBMFT also coincides with large- iV limit of an Sp{N) gen- 
eralization of the SU{2) model.' 

The SBMFT is not meant to provide an accurate quanti- 
tative agreement with 5*17(2) models, nor to replace con- 
trolled numerics. SBMFT should instead be viewed as a 
tool to identify possible phases in competition, their in- 
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stabilities, to provide indicative phase diagrams, and to 
identify the important degrees of freedom that need to be 
incorporated in a description going beyond mean field. 
The most general quadratic Hamiltonian 
would require the use of a larger vector 

(bl^,...,V^^^,b\_^,...,Vj^^^,bi^,...,bNsi) but the 
rotational invariance implies conservation of 

'S'tot = X^i ^If^it + ^li^ii ^^"^ ^lis allowed quadratic 

terms are therefore limited to bl^bjf, bi^b^^^, ^'J^&J^ and 

biibjf. 

With a coordination number 2 and two complex parame- 
ters per bond, it is naively C^'^'. In fact, the modulus of 
self-consistent Aij and Bij cannot exceed a bound, which 
depends on k, see the details in App. B. As for the gauge 
invariance, it allows to fix some parameters to be real, see 
below. 

If we choose the xz plane, we directly obtain a real Ansatz 
from Eqs. 20,21. 



In the most general case, an Ansatz contains MF parame- 
ters for any couple of sites but in practice, we study 
short range interactions. For this example, we limit our- 
selves to first neighbor parameters, but this procedure is 
easily generalized to further neighbors. 
For example, we can deduce Arg(^ij^3jk.4J^i.4im^?mn.4*i) 
from Arg{AijBjkAli) and Avg{AiuiBmnAni)- We can thus 
limit ourselves to fluxes/loops such as each site is en- 
countered exaclty twice. Still this is not enough to have 
independent fluxes, as the fluxes of two loops having a 
common part, with the same operators on the common 
bonds is equal to the flux of the loop encircling both 
loops. For example, we can deduce ATg{AijAj]^AkiAii) from 
Arg(Aj^j*kSki) and Arg(^kiA*iSi*k) (using Bik = SL)- 
If we had chosen a local loop the transformation would be 
simply be a gauge transformation with 8i = tt for all sites 
inside (or outside) the loop and no flux (local or non-local) 
would be modified. 



